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It has recently been shown that spherically symmetric potentials of finite range are uniquely determined by 
the part of their phase shifts at a fixed energy level k 2 > 0. However, numerical experiments show that two 
quite different potentials can produce almost identical phase shifts. It has been guessed by physicists that such 
examples are possible only for "less physical" oscillating and changing sign potentials. In this note it is shown that 
the above guess is incorrect: we give examples of four positive spherically symmetric compactly supported quite 
different potentials having practically identical phase shifts. The note also describes a hybrid stochastic-deterministic 
method for global minimization used for the construction of these potentials. 
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1. Introduction 



be a real- valued potential with compact support. Let R > be a number such 
|x| > R. We also assume that q € L 2 (B R ) , B R = {x : \x\ < R, x G M 3 }. Let S 2 



Let q(x), x El 
that q(x) = for 

be the unit sphere, and a <E S 2 . For a given energy k > the scattering solution ip(x, a) is defined 
as the solution of 



(1.1) AtP + k 2 iP - q{x)iP = , i£l 3 

satisfying the radiation condition 



(1.2) 
where 

(1.3) 

It can be shown, that 



lim 

r — >oo 





dv 


- ikv 


L 




i \ x\—r 


dr 


- V, 


4'Q = 





ds = 0, 
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gifcr /]\ x 

(1.4) ^(x, a) = ipo + ^.(a', a, k) h o I - ) , as r — » oo , — = a' . 

r \r J r 

The function ^4(a/, a, fc) is called the scattering amplitude (0, 

For spherically symmetric scatterers q(x) = q(r) the scattering amplitude satisfies A{a' , a, k) — 
A(a' -a, /c). The converse is established in |Q. Following 0, the scattering amplitude for q = q(r) 
can be written as 

oo I 

(1.5) A(a',a,k) = J2 J2 ^W^MW. 

i=0 m=-I 

where Y/ m are the spherical harmonics, and the bar denotes the complex conjugate. 
The fixed-energy phase shifts — ir < 5i < n (for this k) are defined (fH)]) by the formula: 

(1.6) i4,(*) = ye"'sin(*,). 

Thus, every spherically symmetric potential produces a unique set of phase shifts. The inverse 
problem of finding the spherically symmetric potential q by the set of its phase shifts has long been 
of interest in physics. Details and applications can be found in Newton and Chadan-Sabatier 
(U , where the authors propose a procedure for the identification of the potential q. This procedure 
has been recently examined by A. G. Ramm. (see [1], part 1). 

The following result shows that, theoretically, the inverse problem does have a unique solution 
given a suitable subset of single-energy phase shifts: 

Theorem 1.1 ( ||). Let C be an arbitrary fixed subset of non-negative integers satisfying 

E y = o °- 

ie£,z#o 

Then the data {Si}i e c, corresponding to a potential q(r) 6 L 2 (Bji),R > 0,r = \x\, q = for 
r > R, determine q(r) uniquely. 

Methods for finding piecewise-constant potentials from fixed-energy shifts were developed in 
|lj, Q, ]l5| ]. Numerical results show that, practically, the solution of inverse scattering problem 
with fixed-energy data can be very unstable towards small perturbation of the data. There exist 
quite different potentials which produce practically the same phase shifts (at one energy level). 
The inverse scattering problem with fixed-energy data is ill-posed, see jl| and |fl5|| . After the 
publication of the example in [Q, S. Gutman in [|| has constructed several examples of potentials 
having practically the same phase shifts, and examined the dependency of the phase shifts on the 
energy level k. Based on these observations, a best fit to data method involving several energy 
levels has been proposed and implemented. 

However, the potentials in the above examples change sign. Such potentials can be viewed by 
some physicists as "less physical than one would like to" . It has been conjectured by these physicists 
that it is not possible to construct two different positive potentials which produce practically 
the same phase shifts at a fixed energy and all angular momenta. AGR has discussed this with 
Professor W. Schcid (Giessen, Germany). In this note we construct four positive piecewise-constant 
spherically symmetric potentials having practically identical phase shifts, thus proving that in fact 
such construction is possible contrary to the belief of some physicists. 
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Phase shifts and an algorithm for their computation are discussed in Section 2. A hybrid 
stochastic-deterministic method for finding potentials with practically identical shifts is presented 
in Section 3. Numerical experiments and the constructed potentials are described in Section 4. 

2. Phase Shifts for Piecewise-Constant Potentials 

Below a method from [1] is summarized. Consider a finite set of points = rg < r± < r 2 < 
■ ■ • < rjv = R and a piecewise-constant potential 

(2.1) q(r) — qi, on [r,_i,r,-) for i = 1, ...,N, and q = for r > R. 

Denote nf := k 2 — qi, where i — 1, . . . , N, and k is some fixed positive number. Consider the 
following problem for the radial Schrodinger equation: 

(2.2) ^ + f^- l l+R\ ipi = qipi , hm(2/ + l)!!r- i -VKr) = l, 
which we rewrite as: 

( 2 - 3 ) +\ K i la - I VI = 



on the interval 7"j_i < r < r,. On [r"j_i,rj) one has the following general solution of ( |2.3| ) 

(2.4) W (r) = Aii/^r) + Bini{n t r). 

We assume below that Ki does not vanish for all i. If Kj = for some i then our approach is 
still valid with obvious changes. 

^From the regularity of <pi at zero one gets B\ = 0. Denote Xi = Bi/Ai, then x\ — 0. We 
are looking for the continuously differentiable solution (pi. Thus, the following interface conditions 
hold: 

Aji(niri) + Binifan) = A i+1 ji(K i+1 ri) + B i+1 ni(K l+1 ri), 

(2.5) 

-^[A^Knin) + Bin'^KiVi)] = A l+ ij[(K l+1 ri) + Bj+in^Ki+ir,-). 
The Wronskian W(ji(r),ni(r)) — 1, thus 

Ai+i = nKKi+in^Aijifain) + B t ni(Kir t )} - -^[ni(Ki+iri)[Aifi(Kiri) + B l n[(K i r l % 

(2.6) 

Bl + 1 = 1^3i( K 'i+i r i)[ A ifi( K 'i r i) + B^'^KiTi)] - j'i(K i+ iri)[Aiji(Kiri) + Bini(Kin)]. 
Therefore 

(2.7) M«) = _L(* 

\B i+ iJ Ki+i \a 2 i a 22 J \BiJ 

where the entries of the matrix a 1 can be written explicitly: 

ail = K i+iji( K i r i) n i( K i+i r i) - Kifi(Kin)ni(Ki+in), 

(2.8) 

"21 = «i7'K K * r t)il( K »+l r i) - Ki+lJ(( K i r i)iK K i+l r i); 

a 22 = Kini(Kiri)ji(Ki + iri) - K l+1 ni(Kiri)j[(K i+1 r i ). 
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Thus 

(2.9) 4x+4 2 ^ * 

rv l -L n l w ■ A- 

The phase shift S(k, I) is defined by 

(2-10) w(r) ^LW)l sin(fcr _^ +5(fc)0) r 

where F(k,l) is the Jost function. For r > R 

(2.11) <pi(r) = A N+lJl (kr) + B N+1 m(kr). 



/From ( 2.11 ) and the asymptotics ji (kr) ~ sin(fcr — In/ 2), ni(kr) ~ — cos(£;r — lit/ 2), r — > oo, 
one gets: 

(2.12) tan<5(fc, ?) = - : ^^ ±i = -arjv+i . 



Finally, the phase shifts of the potential g(r) defined in ( |2.l| ) are found from 
(2.13) 6(k, 1) = — arctana;jv+i- 

Let Qo( r ) be a spherically symmetric piecewise-constant potential. Fix an energy level k and 
a sufficiently large N. Let {i$(fc,0}ili be the set of its phase shifts. Let q(r) be another such 
potential, and let {S(k, be the set of its phase shifts. 

The best fit to data function $(q, k) is defined by 



(2.14) *(g,fc) 



sjIii^o-^oi 2 V /2 



For sufficiently large iV such a function is practically the same as the one which would use all 



the shifts in (2.14), since the phase shifts are known to decay rapidly with I, see Jl3j. Our goal is 
to find positive potentials q(r), quite distinct from qo(r), that make the objective function small. 
This is a complex nonlinear minimization problem. An algorithm for its solution is given in the 
next Section. 

3. Global and Local Minimization Methods 

We seek the potentials q(r) in the class of piecewise-constant, spherically symmetric real-valued 
functions. Let the admissible set be 

(3.1) A adm C {(ri,r 2 , ■ ■ ■ ,r M ,qi,q2, ■ ■ ■ ,9a/) : < n < R, qi ow < q m < qugh} , 

where the bounds qi ow and qhigh for the potentials, as well as the bound M on the expected 
number of layers are assumed to be given. 

A configuration (n, r%, . . . , tm, Qi, Q2, • ■ ■ > 9m) corresponds to the potential 

(3.2) q(r) = q m , for r m _ a < r < r m , 1 < m < M , 

where ro = and q(r) = for r > tm = R- 

Note, that the admissible configurations must also satisfy 
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(3.3) n < r 2 < r 3 < ■ ■ ■ < r M ■ 

Given an initial configuration Q 6 A a< i m C R 2M a local minimization method finds a local 
minimum near Q - On the other hand, global minimization methods explore the entire admissible 
set to find a global minimum of the objective function. While the local minimization is usually 
deterministic, the majority of the global methods are probabilistic in their nature. As usual for 
these type of inverse problems the best fit to data function $ has many local minima, see || . In 
this situation it is exceedingly unlikely to get the minima points by chance alone. Thus our special 
interest is for the minimization methods that combine a global search with a local minimization. 
In Q we developed such a method (the Hybrid Stochastic-Deterministic Method), and applied it 
for the identification of small subsurface particles, given a set of surface measurements. The HSD 
method could be classified as a variation of a genetic algorithm with a local search with reduction. 
In jl| two global search algorithms in combination with a special local search method were applied 
to the identification of piecewise-constant scatterers by acoustic type measurements. The global 
algorithms considered are the Deep's Method, and Rinnooy Kan and Timmer's Multilevel Single- 
Linkage Method. The MSLM method has been applied to the identification of piecewise-constant 
spherically symmetric potentials by their phase shifts in || . 

The potentials in this paper have been found using the so-called Reduced Random Search 
Method. 

In a pure Random Search method a batch H of L trial points is generated in A a d m using a 
uniformly distributed random variable. Then a local search is started from each of these L points. 
A local minimum with the smallest value of $ is declared to be the global one. 

A refinement of the Random Search is the Reduced Sample Random Search method. Here 
we use only a certain fixed fraction 7 < 1 of the original batch of L points to proceed with the 
local searches. This reduced sample H ree i of 7L points is chosen to contain the points with the 
smallest 7L values of $ among the original batch. The local searches are started from the points 
in this reduced sample. 

While this method may not be as efficient as the full MSLM algorithm, it has proved to be 
adequate for the problem at hand. 

In our minimization algorithm the Reduced Random Search method is coupled with a deter- 
ministic Local Minimization Method. 

Numerical experience shows that the objective function $ is relatively well behaved in this 
problem. While it contains many local minima and, at some points, $ is not differentiable, stan- 
dard minimization methods work well here. A Newton type method for the minimization of $ is 
described in We have chosen to use a variation of Powell's minimization method which does 
not require the computation of the derivatives of the objective function. Such method needs a 
minimization routine for a one- dimensional minimization of $, which we do using a Bisection or 
a Golden Rule method. See Q| or ||] for a complete description. 

Now we can describe our Basic Local Minimization Method in M 2A/ , which is a modification of 
Powell's minimization method 

Basic Local Minimization Method. 

1. Initialize the set of directions Ui to the basis vectors 

Ui = ei , i = l,2,...,2M. 

2. Save your starting position as Qo- 

3. For i — 1, . . . , 2M move from Qo along the direction m and find the point of minimum Q\. 
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4. Re-index the directions it,, so that (for the new indices) ^(Q[) < ^(<3|) — > ■ ■ ■ > ^K^m) — 

*(Qo). 

5. For i = 1, . . . , 2M move from Qi-i along the direction Ui and find the point of minimum Qi. 

6. Set v = Q 2 m - Qo- 

7. Move from Qo along the direction v and find the minimum. Call it Qo. It replaces Qo from 
step 2. 

8. Repeat the above steps untill a stopping criterion is satisfied. 

Note, that we use the temporary points of minima Q\ only to rearrange the initial directions Ui 
in a different order. 

Still another refinement of the local phase is necessary to produce a successful minimization. 
The admissible set A a dm, see ( 3.1 )-( 3.3 ), belongs to a 2M dimensional minimization space R 2M . 
The dimension 2M of this space is chosen a priori to be larger than 2N, where N is the number 
of layers in the original potential. We have chosen M = 6 in our numerical experiments. However, 
since the sought potential may have fewer than M layers, we found that conducting searches in 
lower-dimensional subspaces of R 2M is essential for the local minimization phase. A variation of 
the following "reduction" procedure has also been found to be necessary in || for the search of 
small subsurface objects, and in Q for the identification of multilayered scatterers. 

If two adjacent layers in a potential have values and Vi and the objective function <E> is not 
changed much when both layers are assigned the same value Vi (or then these two layers can 

be replaced with just one occupying their place. The minimization problem becomes constrained 
to a lower dimensional subspace of R 2M and the local minimization is done in this subspace. 



Reduction Procedure. Let e r be a positive number. 

1. Save your starting configuration Qo — (ri,V2, ■ ■ ■ , Tm, v\, v%, . . . , Vm) € A a d m and the value 
$(Qo)- Let the M + 1-st layer be D M +i = {tm < \x\ < R} and vm+i = 0. 

2. For i = 2, . . . , M+l replace in the layer by Wj. Compute $ at the new configuration 
Qf, and the difference cf = |$(Q ) - $(Qf)|- 

3. For i = 1, . . . , M replace Vi + \ in the layer Di + \ by Vi. Compute $ at the new configuration 
Qf , and the difference cf = |$(Q ) - ^(Q?)i 

4. Find the smallest among the numbers cf and cf. If this number is less than e r Q(Qo), then 
adjust the value of the potential to Vi in the "down" or "up" layer accordingly. Replace the 
two adjacent layers with one occupying their place, and renumber the layers. 

5. Repeat the above steps untill no further reduction in the number of layers is occurring. 

Note, that an application of the Reduction Procedure may or may not result in the actual 
reduction of layers. 

Finally, the entire Local Minimization Method (LMM) consists of the following: 

Local Minimization Method (LMM). 

1. Let your starting configuration be Q = (n, r 2 , ■ ■ ■ , r M ,vi,v 2 , ■ ■ ■ , vm) £ A adm . 

2. Apply the Reduction Procedure to Qo, and obtain a reduced configuration Q 7 containing 
M r layers. 

3. Apply the Basic Minimization Method in ^4 a dmn^ 2M with the starting point Qq, and 
obtain a configuration Q\. 

4. Apply the Reduction Procedure to Q\, and obtain a final reduced configuration Q\. 
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7 



I 


S(k,l) 


I 


5(k,l) 


I 


5(k,l) 





-0.220024E+00 


7 


-0.183339E-02 


14 


-0.204010E-10 


1 


-0.188623E+00 


8 


-0.250850E-03 


15 


-0.766553E-12 


2 


-0.210693E+00 


9 


-0.267137E-04 


16 


-0.253238E-13 


3 


-0.185306E+00 


10 


-0.228367E-05 


17 


-0.741554E-15 


4 


-0.104318E+00 


11 


-0.160476E-06 


18 


-0.193858E-16 


■5 


-0.390310E-01 


12 


-0.944572E-08 


19 


-0.455299E-18 


6 


-0.100159E-01 


13 


-0.472923E-09 


20 


-0.966113E-20 



4. Numerical Results 



Let go be the following potential 



7.2 
4.5 

q (r) = <( 7.2 
4.5 
0.0 



< r < 0.5 
0.5 < r < 1.0 
1.0 < r < 1.5 
1.5 < r < 2.0 
r > 2.0 



Let k = 3. The phase shifts 8(k,l) are computed as in Section 2. They are shown in Table 1. 
Given another piecewise-constant spherically symmetric potential q(r), its ph ase sh ifts 8(k,l) are 
computed in the same way. The objective function &(q ) is formed according to ( 2.14 ) with N = 20. 
This objective function is minimized over A a dm (3.1) with M = 6 and R = 3.0 as described in 
Section 3. A priori bounds for the potential were chosen to be qi ow — 0.0 and quigh = 9.0. The 
Reduced Random Search with the batch size L = 10000 and the reduction factor 7 = 0.01 was 
used for the global minimization. The value e r = 0.1 was used in the Reduction Procedure (see 
Section 3) during the local minimization phase. The initial configurations were generated using a 
random number generator with seeds determined by the system time. The run time was about 30 
minutes on a 333 MHz PC. 

Three potentials with <&((/) < 10 -4 obtained by this procedure are, for example, 



qi(r) 



with = 9.3586605 • 10~ 5 , 



92 0) 



= < 



8 


9991 


< r < 0.4316 




3 


9672 


0.4316 < r < 


8758 


6 


7356 


0.8758 < r < 1 


5718 


4 


3029 


1.5718 < r < 2 


0065 




V 





r > 2.0065 




'6 


4197 


< r < 0.6809 




3 


1509 


0.6809 < r < 


9162 


6 


7464 


0.9162 < r < 1 


2856 


8 


7210 


1.2856 < r < 1 


4314 


4 


5936 


1.4314 < r < 1 


9969 








r > 1.9969 
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with $(q 2 ) = 6.1848208 • 1CT 5 , and 



5.9463 
0.1008 
q 3 (r) = <( 7.9164 
4.6116 
0.0 



< r < 0.8666 
0.8666 < r < 0.9862 
0.9862 < r < 1.4345 
1.4345 < r < 1.9964 
r > 1.9964 



with $(g 3 ) = 3.3089927 • 10~ 5 . 

Figures 1, 2 and 3 show these potentials (solid lines) as well as the original potential qo(r). 
Potentials qi , q2 and q^ are quite different from qo, but their phase shifts are practically identical 
to the phase shifts of go- Since real measurements always contain some noise, such a potential is 
impossible to identify from its phase shifts (at this energy level) by any inverse method. 



5. Conclusions 

Recovery of a spherically symmetric potential from its fixed-energy phase shifts is an important 
physical problem. Recent theoretical results (see Theorem 1.1) assure that such a potential is 
uniquely defined by a sufficiently large subset of its phase shifts at any one fixed energy level. 
However, two different potentials can produce almost identical fixed-energy phase shifts. Such 
examples were obtained in jlj] and ||. These examples give oscillating sign-changing potentials. 
It was suggested by some physicists that such potentials are "less physical". In this note we 
present examples of positive spherically symmetric picccwisc-constant finite-range potentials having 
practically the same phase shifts (at a fixed energy level). The existence of such potentials shows 
that the inverse scattering problem with fixed-energy data can be very ill-posed even if the enrgy 
level is not too small. The solution of the inverse scattering problem becomes much more stable if 
the phase shifts are known at several energy levels, see || for details. It is of interest to study the 
stability of the solution to inverse scattering problem with fixed-energy phase shifts if one assumes 
a priori that the potential is not only positive but also smooth. 

Acknowledgement. AGR thanks Professor W. Scheid for stimulating discussions. 
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Figure 1. Potential qi(r) (solid line), and the original potential qo(r) (dotted 
line); <5>( qi ) = 9.3586605 • 10" 5 . 
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Figure 2. Potential Q2 (r) (solid line), and the original potential qo(r) (dotted 
line); <P(q 2 ) = 6.1848208 • 10~ 5 . 
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Figure 3. Potential qs (r) (solid line), and the original potential qo(r) (dotted 
line); $(g 3 ) = 3.3089927 • lO" 5 . 



